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1. Introduction 

The functioning and development of living organisms - from bacteria to humans - are controlled 
by genetic regulatory networks composed of interactions between DNA, RNA, proteins, and small 
molecules. The size and complexity of these networks make an intuitive understanding of their 
dynamics difficult to attain. In order to predict the behavior of regulatory systems in a systematic 
way, we need modeling and simulation tools with a solid foundation in mathematics and computer 
science. 

In spite of involving a few number of genes, regulatory network usually have complex interaction 
graphs. Many genes are connected to each other and the incoming and outgoing interactions highly 
depend on the gene. This is a major difference with spatially extended systems where a usual 
assumption is that the incoming and outgoing interaction do not depend on the node (translation 
invariance). As in neural networks models, simplifying assumptions are all-to-all interactions 
(globally coupled systems) or random interactions. These assumptions seem not to fit exactly 
with the interaction graphs constructed from biological data analysis. Specific graph models and 
mathematical tools for the analysis of such complex networks of intermediate size are still largely 
missing. 

Genes are segments of DNA involved in the production of proteins. The process in which a protein is 
produced from the information encoded in the genes is called gene expression. During this process 
the genes are transcribed into mRNA by enzymes called RNA polymerase, then the resulting 
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mRNA functions as a template for the synthesis of a protein by another enzyme, the ribosome, in 
a process called translation. The level of gene expression depends on the relative activity of the 
protein synthesis and degradation. In order to adapt the gene expression level to the requirements 
of the living cell, the evolution has selected complicated mechanisms regulating the production and 
degradation of proteins. A simple example of this, is the control of transcription by a repressor 
protein binding to a regulatory site on the DNA, preventing in this way the transcription of the 
genes. Hence, the expression level of a given gene regulates the expression level of another one, 
giving rise to genetic regulatory system structured by networks of interactions between genes, 
proteins, and molecules. 

Let us illustrate the type of interactions used to define our model by means of a simple example: the 
regulation of the expression of the sigma factor a s in Escherichia coli, encoded by the gene rpoS. 
The factor a s takes its name from the fact that it plays an important role in the adaptation to a 
particular stress frequently encountered by bacteria: the depletion of nutriments in the environment, 
which leads to a considerable slowing down of cell growth called the stationary growth phase. The 
protein CRP, a typical repressor-activator, specifically binds the DNA at two sites close to the 
major promoter of rpoS 15 . One of these sites overlaps with the promoter, which implies that 
CRP and RNA polymerase cannot simultaneously bind to the DNA. As a consequence, CRP 
represses the transcription of rpoS. The example illustrates one type of regulation that is quite 
widespread in bacteria: a protein binding the DNA, the regulator, prevents or favors the binding 
of RNA polymerase to the promoter. 

Genetic regulatory networks are usually modeled by systems of coupled differential equations [Hj, 
and more particularly by systems of piecewise affine differential equations like in |Sj. Finite state 
models, better known as logical networks, are also used (see ^H], or HO]). In this paper we 
consider a class of models of regulatory networks which may be situated in the middle of the 
spectrum; they present both discrete and continuous aspects. Our models consist of a network 
of units, whose states are quantified by a continuous real variable. The state of each unit in 
the network evolves according to a contractive transformation chosen from a finite collection of 
possible transformations. Which particular transformation is chosen at a given time step depends 
on the state of the neighboring units. These kind of models are directly inspired by the systems of 
piecewise linear differential equations mentioned before, though they do not correspond to a time 
discretization of those differential equations, but rather to a natural discrete time version of them. 
Such a discrete time model is pertinent when delays in the communications between genes justifies 
discrete time evolution. Besides this, the great advantage of these kind of models is that they are 
suitable for a detailed study from the point of view of the theory of dynamical systems. In this 
framework we may formulate and answer general questions about the qualitative behavior of such 
a system, and eventually relate those results to the experiment. 
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A detailed description of the asymptotic dynamics of our models would require the use of symbolic 
dynamics, as it has been done in jllj for systems related to ours. As a first approximation to 
the complete description of the dynamics we will focus on a global characteristics: the dynamical 
complexity. This is a well-studied notion in the framework of the theory of dynamical systems 
(see ^U], and and it is related to the proliferation of distinguishable temporal behaviors. The 
main motivation of this work is to find explicit relations between the topological structure of the 
regulatory network, and the growth rate of the dynamical complexity. More generally we are 
interested in the possible constraints imposed by the topology of the underlying network and the 
contractive nature of the local dynamics, on the asymptotic dynamical behavior of such kind of 
systems. Hopefully, results on this direction will contribute to the understanding of the observed 
behaviors in the biological systems inspiring our models. The firs steps in this direction has been 
fulfilled by Theorems O O and O below. 

In order to illustrate our results we show some biological motivated examples. Indeed, in sub- 
section w e discuss the self-inhibitor, whose biological prototype is presented in pQ. Then, in 
section we apply our results to particular families of regulatory networks. First we consider the 
family of circuits which models concrete biological examples. For instance, the genetic toggle switch 
studied in |12j corresponds to a circuit on two vertices, while the repressilator in Escherichia coli jH], 
is modeled by a circuit on three vertices. Then, in a biological inspired four vertices network studied 
in j using the results of section |1J we show that the complexity of this particular network may 
correspond to that of a network on two vertices. 

The paper is organised as follows. In the next section we introduce the model, we give the basic 
definitions which we will use in the sequel, and we present some basic facts concerning the relation 
between complexity and asymptotic dynamics. In Section 3 we state and prove the main result 
of this paper, which concerns the complexity of a discrete time regulatory network satisfying a 
(natural but) technical condition which we call coordinatewise injectivity. In Section 4 we explore 
the relationship between the complexity and the structure of the network. Section 5 is devoted to 
some examples, and in the last section we present some concluding remarks. 
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Centre de Physique Theorique de Marseille, and benefited from the financial support of the Min- 
istere delegue a la recherche through the Universite de Toulon et du Var. Both authors thank 
Bastien Fernandez for his remarks, comments, and suggestions throughout this work. 
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To each unit in interaction we associate a label in the finite set {1, 2, . . . , d}. The interaction between 
units is codified by an interaction matrix K E M.dxd([®i !])• The matrix entry Kij specifies the 
strength of the action of the unit i over the unit j. The matrix K is normalized so that, for each j, 
X]f=i Kij = 1- The interaction has two possible modes: activation or inhibition. The interaction 
mode between units is codified by an activation matrix s E Mdxd{{— lj 0, 1}), which is compatible 
with the interaction matrix, i. e., Sij = if and only if Kij = 0. The strength of the interaction 
is weighted by a contraction rate a E [0, 1]; it is directly related to the speed of degradation of the 
state of a noninteracting unit. At each time step t E N, the state of a unit j E {1,2, ... ,d} is 
quantified by a real value x*- E [0, 1], and it changes in time according to the transformation 

d 

x* +1 = ax* + (1 - a)Y J K l , 3 H{s hJ {4 - T M )), 
i=l 

where H : M. — > {0, 1} is the Heaviside function, 

y ' 1 1 otherwise. 

All discontinuities of the transformation are contained in the threshold matrix T E Mdxd([0, 1]), 
which is compatible with the interaction matrix. 

We let the mapping F : [0, l] d — > [0, l] d defined by 

d 

(1) F(x)j := ax* + (1 - a) ^ K^H{s id {y.\ - T id )) 

i=i 

induce a discrete dynamical system on the d-dimensional unit cube, a discrete time regulatory net- 
work. Due to discontinuities, this dynamical system may exhibit a complex behavior characterized 
by the existence of several attractors whose basins of attraction are intermingled in a complicated 
way. Often these attractors consist of a collection of disjoint periodic orbits, but there also exist 
attractors with infinite cardinality. 

For a = 0, the discrete time regulatory network is equivalent to a logical network |13( I19j . For 
a > 0, the model is no longer equivalent to a boolean network, and can be viewed as a discrete 
time analogue of the a system of delay differential equations (see [3] for details). The contraction 
rate a plays the role of (the inverse of) a delay parameter: the smaller a, the stronger the delay is. 
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2.1. The self— inhibitor. In 1 it is studied an autoregulatory system which was experimentally 
implemented to investigate the role of negative feedback loops in the gain of stability in genetic 
networks. This system can be modeled by the one-dimensional regulatory network x i— > F{x) := 
ax + (1 — a)H{T — x), which is semiconjugated to a rotation jlj. As we move the parameters a 
and T, the associated rotation number ranges over an interval, so that the attractor is either a 
finite set, in the case of a rational rotation number, or a Cantor set in the irrational case. Both 
regimes of the self-inhibited unit, periodic and quasiperiodic, share a common feature: the number 
of distinguishable orbits grows at most linearly in time. To be more precise, the orbits {x ,^ 1 = 
F(x°),...,x t = and {y° , y 1 = F(x°), . . . ,x f = F{y t ~ l ), . . .} are distinguishable at 

time t = r if there exists < t < r such that x* and y l lie at opposite sides of the discontinuity, 
i. e., H(T — x t ) = 1 — H(T — y l ). Note that due to the piecewise contractive nature of the dynamics, 
two indistinguishable initial conditions get exponentially close as time goes to infinity. 

The asymptotic dynamics of the self-inhibited unit takes place inside the interval I := [aT, aT+(l — 
a)], where every initial condition get trapped after a finite number of iterations. Inside the attractor 
there are two classes of initial conditions, corresponding to the intervals Iq = [T, aT + (1 — a)] and 
I\ := [aT, T). Since the map x i— > F{x) := ax+(l — a)H(T—x) is injective inside /, and it is an affine 
contraction when restricted to Iq or to I\, then at most one of the intervals F(Iq), F(I\), contains 
the discontinuity. The other one would be completely contained in one of the original intervals 
Iq or I\. Thus, the interval whose image lies inside one of the original intervals contains initial 
conditions whose orbits are indistinguishable at time t = 1, whereas the other interval contains 
initial conditions which may generate up to two distinguishable orbits. Therefore, for t = 1 there 
are at most 3 distinguishable orbits represented by the same number of disjoint intervals. Those 
intervals are precisely the domains of continuity of the mapping x i— > F 2 (x) . We may pursue this 
reasoning and show by induction that for each i 6 N, the maximal cardinality of a set of initial 
conditions in / which generate mutually distinguishable orbits at time t, is at most t + 2. 

For the values of a and T for which the dynamics is semiconjugated to an irrational rotation, the 
number of distinguishable orbits may grow indefinitely. Indeed, for those values, for each n G N 
there exists one interval J d I where F n is a contraction, and such that T € F l (J). This is why 
for these particular choices of a and T there are t + 1 domains of continuity of F t , for each t € N. 
This is precisely the case when the asymptotic dynamics is quasiperiodic, semiconjugated to an 
irrational rotation (see [3] for details). 

In the sequel of this paper we will extend these arguments to the general case, and in this way we 
will obtain upper bounds for the number of distinguishable orbits. 

Let us remark that the system studied in is a generalization of the self-inhibitor we present 
here. 
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2.2. The complexity. In the general <i-dimensional case, given a threshold matrix T G -M<jxd([0, 1]), 
for each coordinate 1 < i < d define the partition V{ of [0, 1] so that the map 

Xj h-> (^(s^^Xj - r i) i)),i?(si i2 (x i - Tj j2 )), • • • ,iJ(si jd (x i - r iid ))) 

is constant restricted to each atom of 7^. We can see that the atoms of Vi are semiclosed intervals 
with endpoints in the set {T i;j : j = 1 . . . , d} U {0, 1}. The map F : [0, l] d -> [0, l] d defined by 
Eq. I^Q) is an affine contraction restricted to the atoms of the Cartesian product V := IIi=i^'i) 
which we call base partition. Indeed, if x,y € [0, l] d belong to the same atom of V, then for each 
1 < hi < d we have H(s itj (xi - T id )) = H(s iij (y i - T itj )), which implies F(x) - F(y) = a(x - y). 

We define recursively the dynamical partitions 

(2) V t+1 := {F' 1 ^) n J : IG?' and J G V}, 

with V 1 := V . It is easy to verify that the r-th iterate of the map F, F T : [0, l] d — > [0, is an 
affine contraction in each atom of V T . It means that two initial conditions in the same atom of V T 
will get closer during the first r time steps. 

The following discussion holds in the slightly more general case when the mapping F : [0, l] d — > 
[0, l) d is a non necessarily affine contraction when restricted to each atom of the base partition V. 
We call this a piecewise contraction. 

The complexity of the system ([0, l] d ,F) is the integer valued function 

C : N -> N such that C{t) := #V l . 

The number of atoms in the dynamical partition V 1 equals the number of distinguishable orbits 
up to time t. It is the aim of this paper to give upper bounds on the growth of this quantity. For 
systems defined by a continuous transformation acting on a compact metric space, if the complexity 
is eventually constant then all orbits are eventually periodic. For the same kind of system, to a fast 
growth of the complexity one can associate some notion of chaos (see [2] and ) . The computation 
of the complexity is part of program whose aim is the complete characterization of the asymptotic 
dynamics of the discrete time regulatory network. Information about the growth of the complexity, 
together with a description of the recurrence properties of individual orbits, will allow us to give a 
detailed picture of the asymptotic behavior of the system. 

2.3. The asymptotic behavior. The limit set A of the system ([0, l] d ,F) is the collection of all 
accumulation points of its orbits, i. e., x G A whenever there exists an initial condition y G [0, 
and a sequence of times t\ < t-i < ■ ■ ■ < t n < ■ ■ ■ , such that hnin^oo F^ tn \y) = x. 

Let us denote 

A := jx G [0, l] d : 3i, k G {1, 2, . . . , d} such that x k = T ijk \ , 
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the discontinuity set. When A n A = 0, then A := Ht^o -^*([P> an d it contains all the orbits 
with infinite past. 

If the limit set has infinite cardinality, then the complexity is a strictly increasing function. Indeed, 
if on the contrary C(r + 1) = C(r) for some r G N, then necessarily to each atom I G P 7 " there 
corresponds a unique atom I such that F(I) C I. The transformation I i-> F(I) := I, acting over 
V T , would generate in this case up to C(r) ultimately periodic orbits. Thus, to each atom IeP T 
we could associate a transient i G N and a period p G N such that F t+kp (I) = F*(I) for each A; G N. 
Since F is a contraction when restricted to each atom of V T , all the orbits with initial condition 
x G I would have one and the same accumulation point in the closure of each one of the atoms 
F*(I), F m (I), • • • , F* +P_1 (I). Therefore, C7(r + 1) = C(r) implies #A < oo. Of course, if A is a 
finite set, then the dynamical complexity has to be a bounded function. 

One would expect bounded dynamical complexity to be a generic situation in the case of piecewise 
contractions. The following suggests this is true. 

Proposition 2.1 (Multiperiodicity) . 7/dist(A, A) > then #A < oo. 

Proof. Let e = dist(A, A)/2 > 0. Compactness allows us to choose a finite cover C = {5(xW,ej) : 
x^ G A, z = 1,2,... N} for A, composed by open balls of radius not greater than e and centered 
at points in A. Since ti < dist(A, A) for each i = 1,2, . . . , N , then F*(.B(xW, e^)) is a ball of 
radius a*ej centered at F*(xW), for each t G N. In this way, for each t G N, the collection 
:= {.B(F*(xW), a*ej) : i = 1,2,..., Af} is an open cover for F*(A), but since A is F-invariant, 
then Ct is an open cover for A. 

Let t be such that a*e is less than or equal to the Lebesgue number of C. Then, for each i = 
1, 2, . . . , N there exists j G {1,2,..., A^} such that 

B(F*(x»),a*e i )cB(xW,e j ). 

For each i = 1, 2, . . . , N choose one such j, and define a map i i— ► r(i) := j from {1,2,..., A^} to 
itself. The successive iterations of this map, i i-> r(i) •-»••••, generate ultimately periodic sequences, 
but the F-invariance of A implies that each i G {1,2, .. . , A^} has to be periodic under r. This 
means that for each i G {1,2,..., A^} there exists p = p(i) G N such that 

F mxip (£(x«, ei )) C B(xW,ei) V m G N. 

Contractivity implies that B{pc^\ ej) H A = {x 1 -^} for each i G {1,2,..., A^}. Furthermore, xW is 
periodic of period (not necessarily minimal) t x Periods are bounded by N xt, so they depend 
on dist(A, A) through the Lebesgue number of C. □ 
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We would like to prove that, in the case of networks of piecewise affine contractions, dist(A, A) > 
is generic. This is true for in the 1-dimensional case [I], but the question remains open for d>2. 

When on the contrary dist(A, A) = 0, the attractor may fail to be invariant under the dynamics. 
In this case, to each point x 6 Awe may associate a "ghost orbits", which would be a real orbit 
by a convenient redefinition of the transformation at the discontinuities. 

Given an orbit O(x) and r G N, let VI := {I G V T : 3 t G N such that x® G I}. This is the 
restriction of the dynamical partition to the orbit. We will say that the orbits O(x) and O(y) are 
distinct if V£ H Vy = for some r G N. Note that distinct orbits are necessarily distinguishable at 
some time. i. e., if O(x) and O(y) are distinct then there exists r G N such that x T and y T belong 
to different atoms of the base partition V. 

We can easily relate the growth of the dynamical complexity to the number of infinite, distinct 
orbits. Indeed, for O(x) infinite, we have #"P* > i + 1 for all t G N. Otherwise there would exist 
t G N such that V 1 = V t+1 and the orbit would be ultimately periodic. Hence, if A contains K 
distinct infinite orbits 0{n i ), 0(x 2 ), . . . , 0(x K ), then there exists a time r G N such that C(t) > 
Y^f =l #P*t > Kit + 1) for each t > r. Therefore, assuming the lowest possible complexity for each 
infinite component of the attractor, a linear behavior of the complexity would imply the finiteness 
of the number of such components. 

The dynamical complexity is directly related to the topological entropy. For a dynamical system 
(X, F), for each r G N and e > 0, a set E C X is said to be (r, e)- separated if for x,y 6 i? there 
exists < t < t such that d(F t (x), F*(y)) > e. The topological entropy of the system (X, F) is 
defined by 

h top := lim(limsuplogA^(r, e)/r), 

e ~ *0 r— >oo 

where N(t, e) = ma~x.{#E : E C X is (r, e)-separated}. 

The topological entropy of a piecewise contraction (X, F) satisfies the inequality 

1 

hop < lim sup - log C (t) . 

t^oo t 

Indeed, for r G N, e > and I G V T fixed, if E C X is (r, e)-separated, then E H I is necessarily 
(1, e)-separated. Therefore N{r,e) < C(r) x N(l,e). 

3. Main Result 

Theorem 3.1 (Polynomial upper bound). Given the interaction and threshold matrices K,T G 
■Mdxd([0, 1]); an d the activation matrix s G M.dxd{{— 1 ; 0, 1}), there exist ao G [0, 1] and c > 1, such 
that for each contraction rate a G [0, ao), the complexity of corresponding discrete time regulatory 
network satisfies C(t) < 1 + c (l + c d t d ) . 
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In the next section we will consider some specific contraints on the structure of the network with 
which we are able to improve this upper bound. We will also show some concrete examples where 
the constant involved in the upper bound is explicitly determined. 

3.1. Prerequisites. For each 1 < j < d let us define the system of affine 1-dimensional contrac- 
tions 



I \i=l / ) 

This system of contractions allows us to compute the coordinates of F. Indeed, for each I G V and 
1 < j < d, there exists an affine contraction / £ T% such that F(pc)j = /(xj) whenever x € I. 

The discrete time regulatory network ([0, l] d ,.F) is said to satisfy coordinatewise injectivity if for 
each 1 < j < d and /, /' £ J-j we have 



To prove the main result of this paper we first show that for networks satisfying this property, C(t) 
grows at most polynomially. In particular, for a network satisfying coordinatewise injectivity, the 
corresponding mapping F is injective. In this case the sets in Q t+1 := F l (V t+l ) are disjoint affine 
images of the atoms in V f , and we can use the equation #"P* = #Q* to compute the complexity. 
This equality holds in the more general situation, when for each t £ N, F 1 maps different atoms in 
V to different sets, not necessarily disjoint. 

In the coordinatewise injective case, the degree of the polynomial bound for the complexity may be 
strictly smaller that the cardinality of the network. In the next section we will see how this degree 
can be related to the topological structure of the network. 

Lemma 3.1. If ([0,l] d , F) satisfies coordinatewise injectivity, then F is injective. 

Proof. For x, y <G [0, l] d with x/y we have two possibilities: either they belong the same atom of 
the base partition, or they belong to different atoms. 

In the first case, since the piecewise constant part of F 



is constant on the atoms of V, it follows that F(x.) — F(y) = a(x — y) ^ 0. 

In the second case there exists 1 < i < d, and different affine contractions /, /' £ T% , such that 
F(x)i = /(xj) while F(y)i = /'(yj. The coordinatewise injectivity ensures that / and /' have 
disjoint images, which implies that -F(x)j / F{y)j, and hence -F(x) — F(y) / as required. □ 



(3) 




/([0, 1]) n /'([0, 1]) = whenever / + f. 
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Proposition 3.1 (Strong contraction implies coordinate-wise injectivity). For each K G Mdxd([®, 1]) 
there exists ao > such that any discrete time regulatory network with interaction matrix K and 
contraction rate a G [0, ao) satisfies coordinatewise injectivity. 



Proof. Fix the threshold and interaction matrices compatible with K. For each j G {1,2, . . . , d}, 
let Tj be the system of affine contractions defined in (j3J). For j G {1, 2, ... , d} fixed, let 

lj : = jr? G [0, 1] : rj := ^eji^j for some e G {0, l} d j . 

Let 5j := min{|r/ — r/| : n,rj' G Ij and n ^ w'}. Finally, let 5 := mini<j<d 5j, which is strictly 
positive. With this define ao := 5/ (1 + 5), so that a < (1 — a)5 for each < a < ao- 

Suppose that a < ao, fix i G {1,2,... , d}, and take /, /' G Ti such that / ^ f. Hence there are 
77, 7/ £ Zj with r) ^ 7/, such that f(x) := ax + (1 — a)r7 and f'(x) = ax + (1 — a)r/. Without loss of 
generality, let us suppose that r? > rj', which implies that rj >rj' + 5. Then we have 

max f'(x) = a + (1 — a)r)' < (1 — a)(n' + 5) < (1 — a)n = min f(x), 
xe[o,i] ' x€[0,l] 

which holds for i € {1, 2, . . . , d} arbitrary, and the proposition follows. □ 

For each 1 < i < d and t G N, let Q* be the collection of intervals 

Q\ ■= { J C [0, 1] : J = IW for some J G Q*}, 

Let us recall that Q t+1 is the image, at time t, of the dynamical partition V t ■ Here and below 
Hj will denote the projection on the i-th coordinate. 

Lemma 3.2. // ([0, l] d , F) satisfies coordinatewise injectivity, then for each 1 < i < d and t G N, 
the intervals in Q* are pairwise disjoint. 



Proof. For each 1 < i < d Qj = V% is a partition, and so the intervals in Qj are pairwise disjoint. 
Suppose that for t = r the intervals in QJ" are partwise disjoint. Let T% be the system of affine 
contractions as defined in equation ©. Then 

Q[ +1 C Ti(QX) V Vi := {/(J) n I : f G ^, J G QJ, and I G 7>i}. 

Since -F satisfies coordinatewise injectivity, and since the intervals in QJ are mutually disjoint, it 
follows that Ti(QJ) V is a collection of pairwise disjoint intervals, implying that the intervals in 
QJ +1 are partwise disjoint as well. □ 
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Given a set of coordinates U C {1,2, ... ,d}, and t G N fixed, a (U,t) -specification is a Tri- 
dimensional interval S G Tligc/ 2i- The bouquet corresponding to the specification S is the collection 
of d-dimensional intervals 

(4) C(S) := {J G Q* : IL;J = ITS V t G 1/'}. 

We will say that J G Q* satisfies specification S if J G C(S), and we will call N(S) := #C(S) the 
degeneracy of S. 

Corollary 3.1. 7/([0, l] d ,F) satisfies coordinatewise injectivity, for each 1 < £,j < d and i G N, 
there exists an interval Jfj G Q\ such that {J G Q* : ly G IljF(J)} C C(Jfj). 

Proof. For each J G Q* we have IljF(J) G J^(Q*) := {/(J) : / G ft and J G Q*}. The previous 
lemma ensures that ^(Q*) consists of partwise disjoint intervals. Thus, for each (i,j) G A, there 
cannot be more that one interval J G Q*, and one affine contraction / £ f;, such that Tjj G 
/(J)- □ 



We denote Jjj the unique interval in Q\ such that if Tjj G IT,F(J) for J G Q l , then ITJ = J\ y 

We say that the d-dimensional interval J G Q* is a predecessor of J' G Q* +1 if J' = F(J) n I for 
some I G T 7 . Reciprocally we say that J' is a successor of J. In our case, since F is injective, every 
J' G Q t+1 has a unique predecessor J G Q*, which we denote by -P(J')- 

Lemma 3.3. Suppose that ([0, F) is coordinatewise injective. Fix U C {1, 2, . . . , d} and t G N. 
For eac/i (U,t + 1) -specification S, there is a unique (U,t) -specification P(S), such that J G C(S) 
imp/ies P(J) G C(P(S)). 

Proof. Fix i £ U and J G Q* +1 satisfying the specification S. Let J' = P(3) and 1 £ V such that 
J = FJ' n I. Then 

(5) ITS = Ui(F(3') n I) = /(IlfJ 7 ) n ikl 

for some / G T%. Since F is coordinatewise injective, then the affine contraction / G T% for which 
(JSJ) holds is uniquely determined from J' G Q*. On the other hand, Lemma 13 . 21 establishes that for 
each 1 < i < d and t G N, the collection Q* consists of partwise disjoint intervals. Hence, for / G T% 
fixed, there is only one interval J' G Q* such that LTjJ' = J' satisfies Equation Therefore ITJ' 
depends only on ITS, and the lemma is proved. □ 

Lemma 3.4. Let ([0, l] d , F) be coordinatewise injective, and let U C {1, 2, . . . , d}. For each t G N, 
and every (U,t + 1) -specification S, we have 

d 

N(S) < N(P(S)) + c x ^ ^ N ( P & x 4i)> 

i0C7 j=l 
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with J* as defined by above, and c := #P — 1. 

Proof. Lemma 13.31 implies that 

C(S) = {P(J') file C(S) : I&V and J' G C(P(S))}. 

If J' € C(P(S)) has more than one successor satisfying S, then necessarily Tj G LLP(J') for some 
i G" £/ and 1 < j < d. Now, each J' G C(P(S)) has no more than #P successors. Therefore 

N(S) < N(P(S)) + c x #{J' G C(P(S)) : T M - G LLP(J') for some i G" U andl < j < d}. 

From Lemma 13.11 if J' G C(P(S)) has more than one successor that satisfies S, then J' satisfies 
P(S) x j| . for some i G" C/ and 1 < j < d. Thus, 

d 

{J' G C(P(S)) : T itj G n,P(J') for some i£ U and 1 < j < d} C J IJ C ( P ( S ) x 4i)' 
and from this we readily obtain the required inequality. □ 

3.2. Proof of Theorem 13.11 Fix a G [0,a ), with a > depending on K, T G A4rf X( i([0, 1]) and 
s G M-dxd({ — 1, 0, 1}) as in the proof of Propositon l3.ll This proposition ensures that the mapping 
F satisfies coordinatewise injectivity. Lemma 13 . 1 1 ensures that F is injective, and so we can replace 
V t by Q* in the computation of the complexity. In order to keep control of the growth of # Q* , we 
use the recurrence Q t+1 = F(Q t ) VP. 

If C(t + 1) > C(t), then the additional contribution in C(t + 1) is due to those ci-dimensional 
intervals in Q* which posses more than one successor. 

Since J G Q* cannot have more than #P successors in Q t+ , one for each possible intersection 
P(J) n I with I G P, it follows that the complexity satisfies the recursion 

(6) C(t + 1) < C(t) + ci x #{3 G Q* : J has more than two successors}, 
where c := #P — 1. 

Now, for J G Q l to have more than one successor it is required that P(J) intersect at least two 
different atoms in P. But if this is the case, then necessarily P(J) contains a discontinuity set 
for P. Hence, if J G Q* has more than one successor, then there exists 1 < i,j < d such that 
Tj G LljP(J). Then, according to © we have 

(7) C(t + 1) < C(t) + c x #{J G Q* : Tj e rLP(J) for some 1 < i, j < d}. 



Using lemma and (J7| we conclude that 

C(t + 1) < C(t) + c x #{J G S* : LLJ = Jjj for some 1 < i,j < d}, 
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and so we can rewrite (JJJ) as 

d d 

(8) C(t + 1) < C(t) + c x £ £ #{J G Q* : n,J = J^}. 

i=l j=i 

Since the cardinalities #{J G Q : LbjJ = j| •} in (jHJ) are the same as the degeneracies N{jj we 
can rewrite (jHJ) as 

d d 

(9) C(i + 1) < C(t) + c x ^ ^ AT^.), 

i=i j=i 

where c := ^P — 1. 

In what follows we will establish polynomial bounds for the growth of the degeneracies N(J), with 
J G Q* arbitrary. For this aim we use a recurrence relation between degeneracies of successive 
specifications. According to Lemma EQ1 for U C {1, 2, . . . , d} and t G N fixed, to each (£7, t + 1)- 
specification S, there corresponds a unique (?7, ^-specification P(S), such that if J G Q* +1 satisfies 
S then its predecessor P(J) satisfies P(S). Then Lemma E3] establishes the inequalities 

_ d 

(10) iV(S) < iV(P(S)) + c x ^ N ( P W x 4i)> 

3=1 

satisfied for each (U, t + l)-specification S, with U C {1,2, . . . , d}. We do not solve this chain of 
inequalities, but instead we solve the chain corresponding to the maximal degeneracies. Indeed, for 
each U C {1, 2, . . . , d} and t G N, let 

n(U,t) := max{A r (S) : S is a (U, ^-specification}. 

Taking the maximum at both sides of the inequality (|1U|) . we obtain the chain of recursions 

(11) n(U, t + 1) < n(U, t) + cx ^(#P; - 1) x n(U U {i}, t) 

taking into account that #{1 < j < d : Tjj > 0} = #Pj — 1. We solved this chain of recursions as 
follows. 

Let us first take U := {1,2,... , c?} \ {i}, for 1 < i < d arbitrary. Taking into account that 
n({l, 2, . . . , d}, t) = 1 for each t G N, from (jllj) we obtain the linear recursion 

n(U, t + 1) < n(U, t) + cx (#P* - 1). 

where c := #P — 1. This recursion has solution n(U,t) < n(U, 1) + c x (#P« — 1) X (t — 1). Since 
n{U, 1) < #Pi < c x (#Vi - 1), we have 

n({7, t) < c x (#7 3 i - 1) x t, V t G N. 
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Suppose that for each U C {1, 2, . . . , d} with #{1 < i < d : i U} = k. Then, 
(12) n(U,t) < c k x t k x Yl(#Vi - 1) V t € N. 

Let [/' C {1, 2, . . . , <i}, with #{1 < i < d : i U'} = k + 1. Substituting this into (fTT|). and taking 
into account the induction hypothesis, we obtain 

n(£/',i + l) < u([f',t) + cx ^(#^-1) xn(U'u{i},t) 

i<£U' 

< n(U',t) + (k + l) c k+1 x t k x J] - 1), 

with solution n{U',T) < n(U f , 1) + (k + 1) c fc+1 x (EI=i **) x Uj^u'(# V j ~ !)> for each r G N - 
Now, taking into account that ^4=1 — r k+1 / (k + 1), we get 

n(£/',T) < c fc+1 x r fc+1 x Yl W?j ~ 1), 

which proves inequality (fT2*)) for each U C {1, 2, . . . , d}. 

Since N(J) < n({j},t) for each t £ N and every J £ Q*, by replacing (fT2|) into ©, we obtain an 
autonomous recursion satisfied by the complexity: 

d 

C(t + 1) < C(i) + c x cfc^ 1 x t d ~ l x - 1). 

i=i 

We solve this recursion taking into account that C(l) = #V = c+ 1 and that n»=i (^^i — 1) < c, 
and finally obtain the upper bound 

C(r) < 1 + ^1 + ^ ^t^ 1 ^ 
< l + c(l + c d r d ) , 

for each r 6 N. □ 



3.3. Generalizations. Theorem 13. II admits two quite direct generalizations as follows. 

Rectangular piecewise contractions. Lemmas from 13.11 to 13.41 can be extended to piecewise 
contractions of the following type: give the base partition V := Yii=i ^i? ^ F : [0, l] d — > [0, l] d be 
such that for each 1 < i < d and I G V, there exists a contractive map fcj : [0, 1] — > [0, 1] such that 

F(x)i = / f)I (xi), VxGl. 

We say that i 7 so defined is coordinatewise injective when for each 1 < i < d, if /ji ^ fo j then 
/;,i([0,l]) n /i,j([0,l]) = 0. Theorem O holds for systems ([0, l] d ,F), where F is a piecewise 
contractions of the type just described, satisfying the coordinatewise injectivity. 
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We say that a piecewise contraction F : [0, l] d — > [0, l] d of the type describe above is a rect- 
angular piecewise contraction compatible with V. We use the adjective rectangular since such a 
transformation, when restricted to an atom of maps the products of intervals into products of 
intervals. 

The bound C(t) < 1 + c(l + c d t d ) holds for any cf-dimensional regulatory network, as far as the 
piecewise injectivity is ensured. Hence, lemma mi implies that given a interaction matrix M, all 
the regulatory networks with interaction matrix M and sufficiently small contraction rate, have 
complexity bounded by the same polynomial 1 1— > 1 + c(l + c d t d ). 

Compatible sequences. Instead of a single rectangular piecewise contraction F : [0, l} d — > [0, l] d 
compatible with a given base partition V, we may consider a sequence SS := {Ft : [0, l] d — > 
[0, 1] of rectangular piecewise contractions compatible with V. Then, we may consider the 
non-autonomous system generated by SS where orbits are given by 

5S (x) := {x* := FtCx'- 1 ), x° = x}. 

The dynamical complexity of the sequence SS is the integer function Css : N — > N such that 

C ss (t) := Wssi where v ss '■= v and Vss '■= { p AQ) n 3 '■ 1 E v ss and J e for t>l. 1£ 
each one of the maps in SS satisfies the coordinatewise injectivity, then lemmas 13. II to HTH hold. In 
this case we use Q^g := F f o ■ • • o F 2 o ^i(^ss) m pl ace °f S*- 
Mutatis mutandi one can prove the analogous of Theorem 13.11 

Theorem 3.2. Let SS := {F t : [0, l] d —> [0,1] }^ be sequence of rectangular piecewise contrac- 
tions compatible with a partition V := ]X=i V a ^ the contractions in SS satisfy the coordinate- 
wise injectivity, then Css{t) < 1 + c(l + c d t d ) with c := #V — 1. 

In what follows a sequence of rectangular piecewise contractions compatible with a given partition 
V, all of them satisfying the coordinatewise injectivity, will be refereed as a V -compatible sequence. 

Remark 3.1. The complexity of a compatible sequence of piecewise contractions with base partition 
V is bounded by the polynomial t *— ► 1 + c(l + c d t d ). Note that the constant c depends only on the 
cardinality of the base partition, indeed, c := j^P — 1. 

4. The underlying network 



A network (directed graph) consist on a finite set V of vertices, and a set of ordered couples 
A C V x V, the arrows. We consider the interacting units to form a network with vertex set V := 
{1,2,..., d} and arrow set A := {(i, j) G V x V : Kij > 0}. We call this the underlying network. As 
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we mentioned above, the main motivation of this study is to understand the relationship between 
the structure of the underlying network and the asymptotic dynamics of the system. 

Let T,K : {1,2,... , cZ} x {1,2, ... ,d} — ► [0, 1] and s : {1, 2, . . . , d} x {1, 2, . . . , d} -» {-1, 0, 1} be 
compatible matrices. Assume that for each 1 < j < d, Ylt=i < 1 an< ^ ^ x a contraction rate 
a G (0, 1). To the sequence of vectors V := {D t G [0, satisfying < D t j + Yli<=v < 1, 

we associate the sequence T>k,t,s,o, '■= {Ft '■ [0, l} d — > [0, l] c! }^ 1 of piecewise affine contractions such 
that 

F t (x)j := a Xj + (1 - a) ( ^ K id H(s iyj (T iyj - x;)) + Aj ) , 

for each i G N and 1 < i < d. Notice that threshold matrix T completely determines the partition 

V := OieV Pi whose atoms are the domains of continuity of the maps Ft G T>T t K,s,a- Taking 
a G (0, 1) sufficiently small so that each transformation F t G T>t : k,s,o, is coordinatewise injective, 
we obtain a "P-compatible sequence. In this way we can define the complexity 

(13) C T ,K, s ,a(t) := max J Cv T ^ s a {t) : V such that < D tJ + K itj < lVt G N j G V \ 

I iev J 

for the parameters K,T,s,a. 

4.1. Skew product. The first result in this direction concerns networks where the dynamics on a 
subnetwork is forced by the dynamics on the complementary subnetwork. Let (V, A) be a network 
such that the vertex set V may be decomposed into two disjoint sets V = Vbase U Vtmndic such that 
A H (Vbundie x ^base) = 0- We say in this case that (A, V) admits base — bundle decomposition. 

When the underlying network (V,A) of the system ([0,l]^ v , F) admits a base-bundle decompo- 
sition V = Vb U Vd, the associate interaction matrix K G A1#y([0, 1]) has the triangular block 
form 

Kb,b Kb,d 

K d4 

where K byb := K\ VbXVb , K b4 := K\ VbXVd and K djd := K\ VdXVd . In this case the projection F(x) 6 := 
{njF(x) : j G V b } depends only on x& := {II^x : i G V b }, thus the projection of ([0, 1]# V ,F) to 
coordinates in Vbase is well defined. 

If the underlying network (V,A) of the system ([0, 1]*^,^) admits a base-bundle decomposition 

V = V b U V d , then the dynamics has the structure of a skew product. Indeed, let 

F b : [0, 1]* V » [0, 1] #V » such that F b (y) := F(y z) 
with z G [0, l]# y<i arbitrarily chosen, and for each y G [0, l]# Vb let 

F y : [0, l]* Vb -» [0, 1]* 1 " 6 such that F y (z) := F(y z). 
Then, for x = x b © x d G [0, we have F(x) = F b (x fe ) © F X6 (xd). 



K = 
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Theorem 4.1. For a regulatory network with parameters T,K,s,a, suppose that the associated 
underlying network (V, A) admits a base-bundle decomposition V = V b UVd- //([0, 1]# , F) satisfies 
coordinatewise injectivity, then 

C(t)<C b (t)xC (Kd>TdjSd>a) (t), 

where t i— > C b (t) is the dynamical complexity of the regulatory network ({0,l]# Vb , F b ), and Kd,Td 
and Sd denote respectively the restrictions of K,T and s to Vd x Vd- 

Proof. Let V b := ILeVi, ^ an d 'Pd '■= Yliev d ^ be respectively the projection of the base partition 
to the coordinates in V b and Vd- Because of the skew product structure of the map, for each i £ N 
and I G V 1 we have intervals G V\ and 1^ G [0, l]# Vd such that I = I& X Id- Since for each 
xj, G [0, l]# Vb , the corresponding bundle transformation F Xb : [0, 1]# — > [0, l]# v depends only on 
the atom J b G V b containing x&, then for each I b G V\ there exists a finite sequence of piecewise 
transformations F u F 2 , . . . , F t : [0, l]* Vd -> [0, l] #yd such 

F'(I t xI (l ) = F 1 oF 1 _ 1 o...oF 1 (I (i ), 

for each multidimensional interval 1^ C [0, l]# Vd . Now, since F is coordinatewise injective, then the 
finite sequence F\ , F2 , . . . , F t can be seen as the prefix of a "P^-compatible sequence in the collection 
T)K b ,T b ,s d ,a £ Qr an a pp r0 p r i a ted sequence of vectors T>. From this it follows that 

N(I b ) := #{IeV t :I = l b x I d , for some I d G [0, l]# Vd } < C Kd , Td , Sd , a (t), 

for each 1^ G V\. From this it follows that 

C(t) = #V l < #Vt X C Kd ,T d ,s d ,a(t) = C b (t) X C Kd ,T d ,s d ,a{tl 

and the theorem is proved. □ 

4.2. Redundant and essential vertices. Theorem VS. II establishes that the complexity of a dis- 
crete time regulatory network, when the contraction rate is sufficiently strong, grows at most 
polynomially. The degree of this polynomial bound is never larger than the cardinality of the ver- 
tex set. In this section we will see how this degree can be related to the structure of the underlying 
network, via the notions of redundant and essential sets of vertices. 

For i, j G V = {1, 2, . . . , d} fixed, for each t G N, and for every J G Q*, let 

(14) Cj(J) = UjC(J) := {ILyJ : J G Q* and HjJ = J}. 

In other words, if the i-th projection of J G Q* is given, i. e., if J satisfies a given ({«'},£)- 
specification, and it is equal to J G Q\, then its j-th projection belongs to the set Cj(J) C Q*. We 
say that i drives j if the cardinality of Cj(J) is uniformly bounded as J ranges over U^Of, i. e., 
if there exists a constant my G N such that #Cj(J) < my for all J G U^ =1 Q t i . 
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A set of vertices U C V is redundant if for each i £ [/ there exists j &V\U such that i drives j in 
the sense stated above. Below we will see why the driving vertex is the redundant one. 

On the other hand, a set W C V is essential if there exists a fixed multiplicity M G N such that 
for each t £ N and S G UieW Qf> iV ( S ) ^ M > L e - ^ C V is essential if the de generacy of 
the (W, ^-specifications is uniformly bounded as t varies. Of course W = V is essential, but the 
interesting case arises when there is a essential proper subset of the vertex set, whose complement 
is redundant. 

Theorem 4.2 (Bounding the degree). Let ([0,l] d ,F) be a discrete time regulatory network satis- 
fying coordinatewise injectivity, and let (V, A) be the underlying network. Suppose that U C V is 
redundant, whereas W :=V\U is essential. Then, there are constants constant c\,C2 > such 
that for each t G N C(t) < 1 + ci (1 + c 2 t q ), with q := #W. 

The proof of this theorem follows the same scheme, with the suitable modifications, as the proof 
of Theorem 13.11 

Proof. Since we are assuming coordinatewise injectivity, Lemma 13 . 1 1 ensures that F is injective, so 
that C(t) = #Q*. As in the proof of Theorem EDI we use the recursion in ©, which can be written 

as 



From here on we use the notion of redundant and essential vertices. Fix U C V redundant, such 
that W := V \ U is essential. Now, since each essential vertex can be driven by several redundant 
ones, let us fix a partition C := {Rj C U : j G W} of U, where to each atom of the partition 
therecorresponds an essential vertex, such that i drives j for each i € Rj. For the rest of the 
essential vertices let Rj = 0. Such partition is not unique since different essential vertices may 
driven by the same redundant one. For C fixed, for each t £ N and j £ W define the collection of 
intervals 



C(t + 1) < C(t) + ci x #{J G Q* : ITJ = J\, } for some (i,j) G A} 
with Jf: as defined by Corollary 13.11 



(15) 




with Ci(Jj k ) as defined by ([Tlji. 

The vertices in Rj are redundant and #{k £V : (i,k) G A} < V% — 1. Therefore 



(16) 
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for each t G N. Here my is the multiplicity bounding the cardinality of Cj(J* •). For a fixed 
partition C and taking into account inequality ©, we obtain 

(17) C(t + 1) < C(i) + Cl x^ ^ iV( J), 
Let us recall that AT(J) := #C(J) = #{J G Q* : I^J = J}. 

As in the proof of Theorem Kill we can find polynomial bounds for the growth of the degeneracies 
N(J). Taking into account the fact that redundant vertices drive essential ones, we may rewrite 
inequality (|lT)|) as 

iY(S)<iY(P(S)) + Cl x E E iV ( P ( S ) xJ )' 

j£W\W> J<£R) 

satisfied for each (W',t + l)-specification S, with W C W. Once again we solve the chain of 
inequalities satisfied by the maximal degeneracies, 

(18) n(W',t + l) <n{V',t) + Cl x ^ rrij n(W' U {j}, t). 

jew\w 

Using exactly the same scheme as in the proof of Theorem 111 II we obtain 

(19) n{W',t) < M c\ t k Yl n i' VteN, 

j£W\W 

with rij : = max(mj, i^Vj). 

Since N(J) < n({j},t) for each t € N and every J € Q*-, substituting inequality <|19|) into (|T7)) . we 
obtain 

cit + 1) < c{t) + C1 e e m tq ' x n ^ 

jew jen* iew\{j} 

< C(t) + M c\ t 9 " 1 E m i IT n i' 

jew iew\{j} 

with q = #W. We solve this recursion by taking into account that rij > rrij and C(l) = c\ + 1, to 
finally obtain the upper bound 

C(r) < 1 + a ( 1 + q M c\ Y[ nj ( E ) 

\ j£W \t=l J J 

< 1 + ci 1 + M c? r 9 n 3 - 
for each r € N. The theorem follows by making C2 = M c\ \\_jew n i- '— ' 
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4.3. Head— independent sets. Now we will exhibit concrete situations where the underlying net- 
work (V, A) admits a non-trivial redundant vertex set with essential complement. In this way 
we obtain concrete examples of discrete time regulatory networks whose dynamical complexity is 
bounded by a polynomial of degree strictly smaller than #V. 

Let us first introduce a characterization of a class of vertices which are essential whenever coordi- 
natewise injectivity holds. 

Fix a network (V, A). A set of vertices U C V := {1, 2, . . . , d} is head-independent if 

1) {j G V : G A} nU = for all i G U, and 

2) {j G V : £A}D{j£V: (i',j) G A} = for all eU with i / *'. 

In other words, U C V is head-independent if 1) each arrow with tail in U has head outside £/, and 
2) different arrows with tail in U have necessarily different heads. We will show that for discrete time 
regulatory network satisfying coordinatewise injectivity, the complement of a head-independent set 
is essential. 

Proposition 4.1. Suppose that the discrete time regulatory network ([0, l] rf ,P) satisfies coordi- 
natewise injectivity. Let (V, A) be the underlying network, let U CV be head-independent, and let 
W := V \ U . Then for each I C V , each t € N, and every (W, t) -specification S we have 

#{J G C(S) : J C 1} < 1. 

Proof. The statement holds for t = 1 since in this case if J G Q} = V is such J C I, then necessarily 
J = I. Now, given a (W, l)-specification S there are two possibilities: 1) either LljS = LTjl for each 
i G W, or 2) there exists i G such that IljSnllil = 0. In the first case #{J G C(S) : J C 1} = 1, 
while in the second #{J G C(S) : J C 1} = 0. 

Let us suppose that #{J G C(R) : J C 1} < 1 for r G N, for each (W, r)-specification R, and for 
every atom I G V. Fix a (W, r + l)-specification S. According to Lemma 13.31 if J G Q T+1 satisfies 
S then P(J) satisfies P(S). Therefore 

{J G C(S) : J C 1} C {F(J')nI: J'GC(P(S))} 

= U VeV {F{3') n I : J' G C(P(S)) and J' C I'}. 

By hypothesis #{P(J') Hi: J' G C(P(S)) and J' C I'} < 1 for each I' G V. Now, for I' G V to 
be such that {P(J') nl : J' G C(P(S)) and J' C I'} + 0, it is needed that IL;I' D n j P(S) for each 
j G VF, which already determines Iljl' for j G W. At the same time, for J' C I' and j G W we 
have IIj(P(J / ) n I) = Sj. Therefore the system of constraints 

(20) ^alLyJ' + (1 - a) K iJ H(s l , j (^ i - T itj ))j fl n 3 -I = n,S, V j G W, 
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must hold for any x£ J'. We consider this system of constraints as a system of equations for the 
unknown binary values {£>jj := H(sij(xi — Tij)) : i £ U, £ A}. Since U is head-independent, 
each one of those equations contains at most one of these unknowns, and each unknown appears 
in at least one of those equations. Coordinatewise injectivity implies that for j £W fixed, the two 
possible values for y^j appearing in 1)20(1 determine two disjoint intervals. Hence, if bij = b is a 
solution for (|2()|). then bij = 1 — 6 is not. Thus, the system of constraints in (12()|) cannot be satisfied 
by two different choices {bij € {0, 1} : i £ U, £ -4.}. Finally, since each projection LTjl' is 

completely determined by the values {i?(sjj(xj — Tjj)) : £ A, Xj £ IljJ'}, then the projections 

rij J', i £ U are determined as well. In this way we fix the remaining projection LTjl' D H J', proving 
that there is at most one atom I' £ V such that #{J' £ C(P(S)) : J' C T, FJ' n I / 0} = 1. □ 

Corollary 4.1 (Complement of a head-independent is essential). Under the same hypotheses as 
in the previous proposition, the set W := V \ U is essential. 

Proof. Indeed, since for each I C V, each t £ N, and every (W, ^-specification S we have #{J € 
C(S) : J C 1} < 1, then iV(S) = £ IeP #{J G C(S) : J C 1} < #P, which proves the claimed 
result. □ 

Remark 4.1. The previous corollary holds if instead of the iterates of a coordinatewise injective 
map we consider a V -compatible sequence of T>x,T,s,a, such that U C V is head-independent with 
respect to the underlying network (V, A) associated to K. 

4.4. 2— loops. In order to have an effective reduction of degree in the polynomial bound for the 
complexity, we need that the underlying network admits redundant sets with essential complement. 
In general, head-independent sets are not redundant, and it seems more difficult to find sufficient 
conditions for vertices to be redundant than sufficient conditions for vertex sets to be essential. 
Below we will show a case where the structure of a two-vertices subnetwork is such that one of its 
vertices drives the other. 

We say that the interaction matrix K £ Mdxd([0, 1]) is non- degenerated if for each j £ {1,2, ... , d}, 
the map 

d 

je £ {0, l} d : e» = for all i such that K itj = oj 3 e i-> ^ ejJQj £ [0, 1] 

i=i 

is one-to-one. 

Fix a network (V,A). A 2-loop in (V, A) is a couple {i,j} C V such that £ A and i is the 
only tail for an arrow whose head is j, and ( j, i) £ A and i is the only head for an arrow whose tail 
is j. In this case we say that j is the isolated end of the 2-loop. 
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Proposition 4.2. Suppose that ([0,l] d ,F) satisfies coordinatewise injectivity, and suppose that 
C V is a 2-loop in the underlying network (V, A) whose isolated end is the vertex j. If in 
addition the interaction matrix K is non-degenerated, then i drives j. 

Proof. We will prove that, for each t E N, each J E Q*, and every / E Vj, there exists at most one 
interval J' E Q* such that J' C I, and C(J x J') / 0. In this way we prove that i drives j with 
multiplicity rriij < jf-Vj = 2. 

For t = 1 the statement obviously holds. Suppose it holds for t = r. Fix J E QJ +1 and / C Vj. 
Then, thanks to the coordinatewise injectivity, there exists a unique contraction / E Ti and a 
unique interval P(J) E QJ such that J = f{P(J)) n I 7 . Here I' is the atom Vi which contains J . 

The affine contraction / : [0, 1] — > [0, 1] has the form /(a;) = ax + (1 — a) X^fc=i e kKk,i, for some 
e E {0,l} rf . Since the interaction matrix is non-degenerated, then the values of E {0,1}, for 
Kj- i > 0, are uniquely determined by /. Therefore, given J E Q[ +1 , if J' = P(J) for some J E C(J), 
not only we have J' E C(P(J)) as established in Lemma 13.31 but also H(ski(^-k ~ ^fei)) = e fc> f° r 
each x E J'. 

From here on we use the fact that {i,j} is a 2-loop with isolated end j, i. e., (a) Kj k = for 
k ^ i, and (b) iffej = for k ^ i. Condition (a) implies that the value ej E {0, 1} determines the 
atom I" E Vj containing HjJ' for each J' E C(P(J)). Here we apply the induction hypothesis. 
Given P(J) E QJ and 7" E Pj, there is at most one interval in QJ, which we denote J', such that 
J' C I" and C(P(J) x J') ji 0. By condition (b), the image of each J' E C(P(J) x J') is the 
interval UjF(J') = aJ' + (1 — a)iJ(sjj(xj — Tjj)), which depends only on J and -P(J) 9 Xj. Thus, 
the interval Hf-FfJ 7 ) fl J E QJ +1 is uniquely determined from J E Q^ +1 and I £ Vj. □ 

Remark 4.2. Once again, proposition \4-°A holds for compatible sequence of maps with common 
base partition, all of the satisfying piecewise injectivity, and having the same underlying network. 

5. Examples 

5.1. Circuits. The circuit on d vertices has vertex set V = {1, 2, . . . , d} and arrow set A C V x V 
satisfying 

(i, j) E A<$ j = i + 1, or i = d and j = 1. 

Some numerical experiments suggest that in large random networks, and for a large class of initial 
conditions, the state of most of the units in the network converges to a fixed value, whereas in 
a small subnetwork, an oscillatory regime takes place. Very often this subnetwork is a circuit of 
low dimension. The circuit on d vertices may be negative or positive according to the value of the 
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product a := s^i X Yii=i s «,«+ii defined by the activation matrix. In discrete time regulatory net- 
works defined in positive circuits, numerical experiments have shown mutistability, i. e., coexistence 
of multiple attractive fixed points. On the contrary, for the systems defined on negative circuits, 
beside some degenerated cases, the dynamics does not admit fixed points. A rigorous study of this 
numerically observed phenomena is in progress. Concerning the complexity of this kind of systems 
we can already state the following. 




Figure 1. 



A discrete time regulatory network whose underlying network is a circuit on d vertices, and whose 
contraction rate is smaller than 1/2, satisfies C(t) = 0{t d ). This follows from Theorem 13 .![ taking 
into account that these systems satisfy coordinatewise injectivity. Indeed, since the interaction 
matrix is normalized, i. e., X^=i-^i,j = 1j then K{j = 1 for each € A. Using the same 

notation as in Proposition 13.11 we have that in this case It = {0,1} for each i € {1,2, ... ,d}, 
implying that ao = 1/2. 

The genetic toggle switch described in |12| is a concrete experimental example of a positive 2- 
circuit, made of two mutually inhibiting interacting units. In this case each unit is composed by 
a couple promoter-repressor, such that each promoter is inhibited by the repressor transcribed by 
the opposite promoter. This system is traditionally modeled by the system of ordinary differential 
equations 

dr\ ol\ 

IT ~ i + r A " ri 

dr 2 __ a 2 

where rj represents the concentration of repressor i, ai is the effective rate of synthesis of repressor 
i, and (3i is the cooperativity of repression of promoter i. The corresponding network of piecewise 
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contractions has the form 



n(t + 1) = on(t) + (l - o)fr(r lj2 - r 2 (t)) 

r 2 (t + 1) = ar 2 (t) + (1 - a)iI(T 2il - n(i)), 



where the rate of synthesis of the repressors is controlled by the parameter a, and the cooperativity 
of repression of promoters is related to both a and the thresholds TJj. The dynamical behavior 
of the positive circuit in two vertices is completely understood (see [3]), and it can be either 
bistable (the same as the continuous model) or conjugated to a rotation. In all cases the dynamical 
complexity asymptotically bounded by t h- »■ t + 1. 



For the general circuit on 2 vertices, if the contraction rate is smaller that 1/2, each one of the 
2 vertices is at the same time redundant and essential. Indeed, in this case the coordinatewise 
injectivity holds, each vertex constitutes a head-independent set, and at the same time it is the 
isolated end of a 2-loop. Corollary 14.11 and Proposition 14.21 apply, implying that each vertex 
is redundant and essential at the same time. Hence, a system of two interacting units whose 
underlying network is a circuit on 2 vertices, and whose contraction rate is smaller than 1/2, 
necessarily has linear complexity. Furthermore, for the negative circuit, since the image of rectangle 
in Q T cannot intersect the double discontinuity (xi = Xi )2 ,x 2 = T 2) i), then it cannot have more 
than one successor. Taking into account this constraint in the computations we performed above, 
we deduce that in this case C(t) < 2t + 2 for each f £ N. It has been proved [3] that the 
negative circuit admits an arbitrary large number of periodic orbits, as one moves on the line 
{(T, a) : T12 = T%,\ = 1/2, a € [1/2,1)}. There the number of admissible periodic orbits grow as 
one increases the contraction rate. Because of this, the complexity increases as a grows, though our 
numerical experiments indicate that C(t) can always be bounded by a quadratic polynomial, as it 
is illustrated in Figure 2. We conjecture that the complexity of the negative circuit on 2 vertices 
can always be bounded by a linear function C(t) < K X (t + 1), where the constant factor depends 
on the choice of parameters (T,a), and it diverges as we increase a. 
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Figure 2. 

Log-log plot of the complexity of the negative circuit on 2 vertices 
for Ti ; 2 = T<2 \ = 0.5 and a = 0.93. The straight line is the log-log 
plot of the quadratic upper bound t — » 4i 2 . 

For a circuit on d vertices any independent set of vertices (i. e., no two of them are adjacent to the 
same arrow) is head-independent. In particular {2,4,... ,2[<i/2j} is a maximal independent set, 
therefore, according to Corollary 14. II its complement {1, 3, . . . , 2\_{d — l)/2j + 1} is essential. If we 
think of the collection Q T as a finite time version of the attractor, Proposition 14. II establishes that 
this coarse grained attractor has dimension k = \_d/2\ . Indeed, in this case each (f-dimensional 
interval of J G Q T has at most k independent one-dimensional projections, whereas d—k projections 
are determined by the first, and by the atom of V containing J. In this way, the coarse grained 
version of the attractor lies inside the graph of a #"P-valued function defined on [0, l] fc , and taking 
values on [0, l] d ~ k . Note that this property holds uniformly on r, suggesting that the attractor 
itself may be embedded in the graph of a multivalued function. This of course has yet to be proved. 

Let us remark that in j^j it is presented an biological example of a negative circuit on 3 vertices, 
the so-called repressiliator. They report the existence of oscillatory behavior in time, which would 
correspond to periodic solutions of the corresponding model, though the period and amplitude 
exhibit significant variability. This may correspond either to stochastic effects, or may reflect 
intrinsic complex behavior. 

5.2. Networks with loops. With respect to the effective reduction of the polynomial upper 
bound, the example of the circuit on 2 vertices with contraction rate smaller than 1/2 may be 
generalized as follows. Fix a network (V,A) and 2-loops {i,j} and {i',j'}, with isolated ends j 
and j' respectively. We say that and {i',j'} are disjoint" if {k G V : (i, k) G .4} PI {k G V : 

(i', k) G A} = 0, i. e., the non-isolated ends have disjoint heads. 
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Let (F, [0, 1] ) be a network with non-degenerated interaction matrix (as denned in the previous 
section), and satisfying the coordinatewise injectivity. If the underlying network (V,A) has k := d 
partwise disjoint 2-loops, then C(t) = 0{t d ~ k ). Indeed, let {{i n ,jn} ■ n = 1,2, . . . ,k} be the 
set of mutually disjoint 2-loops, and let j n is the isolated end of {i n ,jn} for each n = 1, 2, . . . , k. 
According to Proposition 14,21 the set U := {i n : n = 1, 2, . . . , k} is redundant. Indeed, each j n 
drives i n for each n = 1, 2, . . . , k. On the other hand, since the 2-loops are partiwise disjoint, the 
set U head-independent too, and so by Corollary 14.11 the complement W := V \ U is essential, 



In Figure 3 we represent a network containing 3 disjoint 2-loops with isolated ends a, b and c. 
According to Theorem 13.11 and Corollary 14.21 if the contraction rate is sufficiently small and the 
interaction matrix is non-degenerated, then the complexity of the system is bounded by a cubic 
polynomial. 

In a discrete time regulatory network of the type considered in this paper was used to model the 
regulation of the expression of the tumor suppressor gene p53, whose product acts as an inhibitor 
of uncontrolled cell growth. The underlying network of this model is shown in Figure 4. The node 
m represents the expression level of the mdm2 gene, whose product inhibits the expression of p53. 
Nodes b and c represent blood-vessel formation and cell proliferation respectively. Though b and c 
are not expression levels of a particular gene, they represent relevant quantities directly implicated 
in cancer decease. 



therefore C(t) = 0{t d - k ). 





Figure 3. 
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Figure 4. 

Regulatory network for the p53 gene expression. The arrow indicates 
positive interaction, while H indicates an inhibitory interaction. 

In this network the couple p53-m form a 2-loop with isolated node m, and the vertices m and c are 
head-independent. If we assume coordinatewise injectivity, according to proposition ^. 11 nodes p53 
and c capture the whole complexity of the network. On the other hand, since node m is redundant, 
the complexity of the network is bounded by a polynomial of degree 3. Note that in this case the 
network admits the base-bundle decomposition {p53,m} U {&, c}, and both base (V& := {p53,m}) 
and bundle (Va '■= {b, c}) define circuit on two vertices. Therefore, if coordinatewise injectivity 
holds, Theorem 14.11 implies that the complexity of this four-vertices network is bounded by a 
quadratic polynomial. 

6. Final Remarks 

In this work we formulate and show the first results towards the solution of the following theoretical 
problem: determine to what extent the topological structure of the underlying network determines 
the behavior of a discrete time regulatory network. By solving this problem we intend to acquire 
a better understanding of the behavior observed in real life regulatory networks. The next step 
in our program consist in extending Theorem 13.11 to the case of arbitrary contraction rates, when 
coordinatewise injectivity does not hold. Some preliminary numerical experiments support lead us 
to conjecture that the complexity is always bounded by a polynomial; furthermore, this experiments 
suggest that the asymptotic dynamics consist of a finite number of periodic or quasiperiodic orbits. 
If this is true, the complexity should be always bounded by a polynomial of degree one. 

Though in general it is not possible to deduce relevant dynamical properties only from the com- 
plexity, we expect that in our case the growth of the complexity could give information about 
topological and recurrence properties of the attractor. For instances, the growth rate of the com- 
plexity could be related to the number of transitive components of the attractor. Anyway, in order 
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to have a complete description of the qualitative behavior of our models we need to consider other 
characteristics besides the dynamical complexity. 

Theorems 14.21 and 14.11 we have specific instances of constraints imposed on the behaviors of the 
system by the structure of the underlying network. In this direction, our next goal is to study 
how the behavior of a network is related to the possible dynamical behaviors of its subnetworks. 
Results in this direction would allow us to deal with questions such as the role of self-regulation in 
the stabilization of the circuits, and in general, to progress in the understanding of the engineering 
behind the design of real life functional networks. 
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